Analytical solutions to zeroth-order dispersion relations of a cylindrical metallic nanowire 
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Zeroth-order complex dispersion relations of a cylindrical metallic nanowire have been solved out analyti- 
cally with approximate methods. The analytical solutions are valid for the sections of the dispersion relations 
whose frequencies are close to the Surface Plasmon frequency. The back bending of the Surface Plasmon- 
Polaritons(SPPs) can be well described by the analytical solutions, confirming that the back bending is orig- 
inated from the metal Ohmic loss. The utility of the back bending point in the dispersion relation for the 
measurement of the metallic Ohmic loss has also been suggested. 
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I. INTRODUCTION 

Surface Plasmon-polaritons(SPPs) induced by the coupling 
between the free-electron plasma and electromagnetic field at 
the metal-dielectric interfaces are known for the wide appli- 
cation in nano optics. 15 Due to the spatial localization at the 
metal-dielectric interfaces, the SPPs can be guided to manip- 
ulate light in nanoscaled photonic circuitry. 2-4 It has also been 
found other effects of SPPs, including field enhancement 6-8 , 
and influence on lifetime of emitters close to the SPPs . 9 ' 10 
In various plasmonic structures, cylindrical metallic nanowire 
— — as a symmetric plasmonic structure has been widely stud- 
ied for the applications of the optical waveguide and the scan- 
ning near-field optical microscopy (SNOM). In order to un- 
derstand the physical properties of SPPs in those structures, 
it is important to get the dispersion relations of the frequency 
and wave vector of the SPPs, which determine the basic prop- 
erties of the SPPs. Principally, the dispersion relations can be 
obtained by solving the Maxwell equations with the boundary 
conditions of the plasmonic structures imposed. For the cylin- 
drical nanowire, the SPPs dispersion relations can be numer- 
ically calculated from the dispersion equationjiiii The SPPs 
dispersion relations of plasmonic structures are complex when 
the metal Ohmic losses are introduced] 15,17,18 There exist two 
types of solutions to the complex dispersion relations. One 
solution specifies a complex frequency as a function of a real 
wave vector, noted as complex -co for convenience. The other 
one specifies a complex wave vector as a function of a real 
frequency, noted as complex - k. Even though the two types 
of solutions are obtained from one same structure, they are 
rather different. For the sections of the dispersion relations 
with frequencies close to the surface Plasmon frequency, the 
complex-k solution has a back bending, which is absent in the 
complex - oj solution. The latter is an asymptotic curve.^^ 
The back bending of the dispersion relation has some 
unique properties. For example, It has been calculated that 
the minimum energy velocity is exactly at the back bending 
point.i^— However, the utility of the bending back behavior 
of the SPPs to control light has not been fully explored yet. 
Furthermore, the local density of states(LDOS) of SPPs is a 
key parameter influencing the lifetime of an emitter close to 
the SPPs . 9 >— According to the definition of the LDOS that is 
proportional to 1 /|V^<y(A:)| , the LDOS close to the back bend- 
ing point in the complex - k dispersion has a small and broad 



peak, while the LDOS in complex - co dispersion has a large 
and narrow peak due to the asymptotic curve. Thus, in or- 
der to compare the physical results between the two types of 
solutions, it is necessary to understand the behaviors of the 
asymptotic curve and the back bending of the SPPs. 

It has been suggested that the complex-k solution describes 
the SPPs mode decaying spatially while the complex - co so- 
lution is for the SPPs decaying in time rather than in space 
i 17 i 18 i 21 The discrepancy of the two solutions has been consid- 
ered to be originated from the metal Ohmic loss . 22 ~ 25 Such 
conclusion is originally drawn from the study of SPPs on pla- 
nar metallic-dielectric surfaces. When a perfect metal with- 
out damping is considered in this planar case, the two types 
of solutions of the SPPs dispersion relations overlap with the 
same asymptotic behavior. When the metal Ohmic loss is 
introduced into the dielectric response of the metal, the two 
types of solutions are then different. Before such conclusion 
is extended to the case of the cylindrical metallic nanowire, 
it is necessary to exhibit how the metal Ohmic loss leads to 
the back bending of the dispersion relations of the metallic 
nanowire. In this paper, we derive approximate analytical 
solutions to the complex dispersion relations of a cylindrical 
metal nanowire for the description of the back bending as well 
as the asymptotic behavior of the SPPs. Based on the analyti- 
cal expressions, the role of metal Ohmic loss in the determina- 
tion of the back bending is discussed. The utility of the back 
bending for the measurement of the metal Ohmic loss has also 
been suggested. 



II. DISPERSION EQUATION 

The metal nanowire is modeled to have the shape of a cir- 
cular cylinder with the radius r and the infinite length in a 
medium. The electromagnetic field in this model can be ex- 
panded with cylindrical harmonics. By imposing the elec- 
tromagnetic boundary conditions at the metal-dielectric inter- 
face, we obtain the following transcendental equation: 
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which repeats the reported results. 11 '— Here, is the first 
kind of Hankel function with the order of integer n and J„ 
is the Bessel function of nth order. The Hankel and Bessel 
functions with denote represent the first order differentiation. 
kj(j=o,i) is the wave vector with the value of kj = yejw/c, 
where £ 7 is the dielectric function and the subscript j labels 
the quantities outside the nanowire (j = 0) or inside it (j = 1). 
c is the speed of light in vacuum. The dielectric function of 
the metal can be expressed as 



Ei(w) = £«,[! 



a)(a> + z't) 



], 
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where co p is the bulk-plasmon frequency and t is the bulk elec- 
tron relaxation rate, 15 which reflects the metal Ohmic loss. 800 
in equation (2) is a constant for the general description of the 
dielectric function of the metal. k z is the component of the 
SPPs wave vector along the cylinder axial and the radial com- 
ponents of the wave vectors are defined as k r j = Jkl — k%. We 

note that the value of k r o should be chosen to guarantee the 
imaginary part of k r Q to be positive since the light intensity 
should be decaying away from the metal cylinder. The dis- 
persion relation between k z and a> then can be obtained from 
the equation (1). However, it is impossible to get an analytical 
solution to the equation for the whole dispersion relation. We 
show in the following that the sections of the dispersion rela- 
tions whose frequencies are close to the Surface Plasmon fre- 
quency can be solved out from equation (1) analytically with 
approximate methods. For clarity in the following derivation, 
we renormalize the a> and r by aj p . Wave vector components 
of kj(j = o,i) , k z and k r j(j=o,i) are renormalized by a> p /c and r is 
renormalized by c/co p . 



III. ANALYTICAL SOLUTIONS 

For the sections of the dispersion relations with frequencies 



close to the Surface Plasmon frequency u sp = J ^| , the 

SPPs modes of the metal nanowires are nonradiative and the 
real parts of k, noted as Re[k z ] at those sections have larger 
values than the imaginary parts Im[k z ]. The values of Re[k z ] 
play more important role in the Bessel and Hankel functions 
than the values of Im[k z ], which makes the Bessel functions 
and Hankel functions in the equation (1) behave as modi- 
fied Bessel functions. In order to verify this point, we have 
conducted numerical calculations on the Bessel functions and 
Hankel functions with the frequency close to the Surface Plas- 
mon frequency, and found that the Bessel function approaches 
to infinite while the Hankel function goes to be zero. Thus, it 
is reasonable to replace the Bessel and Hankel functions in 
equation (1) by the following approximate expressions: 
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with t replaced by -ik r jr in the derivation. In order to get 
the analytical solutions to equation (1), fcr/CMU) should be ex- 
panded by k z as 



k 2 



2k 2 



2L 



(3c) 



In this paper we only discuss the zeroth-order analytical solu- 
tions with n — 0. Substituting equation (3) into the equation 
(1), we obtain the following equation 
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There exist four solutions to the equation (4). However, only 
one solution is correct. Especially for the case of r — > 00, the 

9 k 2 k 2 

equation (4) is then solved out to be kz = jt-jt, which is the 
right dispersion relation of SPPs on a planar metal surface. 



A. Solutions for e = 1 and = 1 

For clarity to show our derivation, we set the dielectric con- 
stants as 80 = E M = 1 . Then, the Surface Plasmon frequency 
is u> sp = 1/ v2~. Firstly, we consider the complex - k solution. 
We rewrite the equation (4) expressed by k z , a>, t and r as 



[(4a>*-2) + (*r)(- 



)+4(iT) z W z + [- + — ]kf 



a> ~ r ru 

+[6w 2 (l - w 2 ) - 1 + 6w(/t)(1 - 2lo 2 ) - 6co 2 (iT) 2 ]kl 

+(2u> 2 - l)(a> 2 - l)co 2 + (/t)(4w 2 - 3)w 3 + 2(rr)V 
= 0. (5) 

In the complex - k solution, there exists a back bending of the 
SPPs where the frequency 10 is close to the Surface Plasmon 
frequency cj sp — 1 / y2. Comparing to the terms with the same 
order of it, we find that the last three terms without k, can be 
neglected since the Re[k z ] is large enough at the back bending 
point. Thus, equation (5) can be further simplified as 



7 8ar — 2 *, 1 (if) 
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+[6w 2 (l - to 2 ) - 1 + 6w(/t)(1 - 2u 2 ) - 6co 2 (iT) 2 ] 

= 0. (6) 

Then, the complex k z can be solved out from the equation (6). 
The approximate methods in above derivation are valid for 
larger Re[k z ] than Im[k z ], which has been confirmed at the 
sections of the dispersion relations by our numerical calcu- 
lations. And in the derivation, the approximate methods are 
conducted regardless of Im[k z ]. Thus, the correct solution to 
the equation (6) for the Re[k z ] is 



(3b) 



Re[k z ] = Re[ 
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with 

k A = 4to 2 - 2 + ) + 4(;t) 2 , 

to 

k c = 6w 2 (l - to 2 ) - 1 + 6w(/t)(1 - 2w 2 ) - 6w 2 (/r) 2 . 
Especially for the case of t — 0, /?e[fe»] is then equal to 



-i- Ji- 4(4 W 2 -2)[6 W 2 (l-w 2 )-l] 

= " 2 „ ■ (8) 

2{4to A - 2) 

Now we consider the complex - to solution. We rewrite the 
equation (5) expanded with the frequency to as the following 
equation: 

2a) 1 + 4(/r)w 6 - [6k 2 - 2(/r) 2 + 3]to 5 - [12(/r)/fc 2 + 3(/t)]« 4 
+[4kt - 6(ir) 2 k 2 + 6k 2 + + [S(ir)kt + 6(ir)k 2 ]to 2 



-[-4(h) 2 k 4 z + 2k\ + k\ - k 3 Jr]o) + (/r)fc?/r - 2(h)k 4 z 
--0. 



(9) 



Similarly, comparing to the terms with the same order of ir, 
we simplify equation (9) by neglecting the terms with small 
values to get 

to 3 + 2(iT)u 2 + [l/(4rjy - 1/2] w + (rr)[l/(4rJfc z ) - 1/2] = 0. 

(10) 

The imaginary part of to has been calculated to be small [15]. 
Thus, the correct approximate solution to equation (10) is 



Re[to] = Re[ 



(ii) 



with 
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B. General solutions 

For the general case with the dielectric constants eo and 8.x, 
introduced into the wave vectors, analytical solutions to Equa- 
tion (4) can be obtained similarly by neglecting the terms with 
small values. The approximate methods are still valid for the 
sections of the dispersion relations with frequencies close to 

the Surface Plasmon frequency u sp = J ^+ z ■ Here, we only 

show the results. For the complex — k solution, the Re[k z ] is 



Re[k z ] = Re[ 
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with 
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For the complex - to solution, the Re[co] has the same forma- 
tion to the equation (1 1) as 



Re[to] = Reb 



( i )J+( | )3+ ^_i + ^ +( | )3] , 

(13) 



with 
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IV. DISCUSSION 

The solutions of equation (7), (11), (12) and (13) are the 
main results of this paper. In order to verify these results, we 
perform direct numerical calculations to solve the equation ( 1 ) 
for the comparison (Figure (1)). The comparison is conducted 
on an Aluminum nanowire with the radius of 100«m, which is 
corresponding to the renormalized parameter r = 1.90415. 16 
The renormalized r takes the value of r = 0.039062 and di- 
electric constants of Eo = £«, = 1 are for the comparison. 
Figure 1(a) is for the comparison of the complex - k solution 
while figure 1(b) is for the complex-to solution. The scattered 
data in the figure are the exact solutions to the equation (1) by 
the direct numerical calculations, while the solid line in fig- 
ure 1(a) (figure 1(b)) is obtained from equation (7) (equation 
(11)). It shows that the Re[k z ] at the section of the complex— k 
dispersion relation close to the back bending point in figure 
1(a) can be well fitted by the equation (7) and the asymptotic 
curve of Re[oj] in figure 1(b) can be described by the equa- 
tion (11) accurately. We have also conducted the comparisons 
with various r and r for the Aluminum nanowires and get the 
same result of figure 1 . To show the result, we only summa- 
rize the comparisons of complex -k solutions in figure 2 since 
the comparisons of complex - to solutions are tedious with the 
same asymptotic curve. We plot k z o, the value of Re[k z ] at the 
back bending point indicated in figure 1 (a), as functions of the 
r and r in figure 2. The scattered data are the exact solutions 
to the equation (1) by direct numerical calculations while the 
solid lines are the analytical results from the equation (7). In 
figure 2(a), the t is fixed to be r = 0.039062, while in figure 
2(b) the r is fixed to be r = 0.2. The coincidence between 
the scattered data and the solid lines in figure 2 confirms that 
our analytical solution of equation (7) can be used to accu- 
rately describe the section of SPPs dispersion close to the back 
bending point. 

Figure 2 also shows that k z o is sensitive to the radius and the 
metal loss of the nanowire, especially for the nanowire with 
radius less than 50nm. For the SPPs application, the metallic 
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FIG. 1: dispersion relations of a metallic nanowire with the shaper 
of a circular cylinder, (a) complex - k solution to the relation of 
Re[k z ] and frequency a>. (b) complex - co solution to the relation 
of Re[u>] and k z . Scattered data are the results from equation (1) by 
numerical calculation. Solid lines in Figurel (a) and (b) are obtained 
from equations (7) and (11) respectively. For the calculation, the 
renormalized parameters r and t take the values of r = 1.90415 and 
t = 0.039062. Dielectric constants take the value of e = 1 and 
e„ = 1. 



nanowires are fabricated normally with the radius having the 
order of few nanometers. The size of such nanowires can be 
identified by the Atomic Force Microscopy (AFM) and Trans- 
mission Electron Microscopy (TEM) techniques. Thus, we 
suggest that equation (7) can be used to fit the metal loss of 
the nanowire with the radius given by AFM or TEM, if the 
back bending of the dispersion relation can be experimentally 
obtained. In this way, the difference of the metal loss between 
the bulk material and the nanowire can be distinguished. It has 
been reported that the attenuated-total-reflection techniques 
can be used to measure the SPPs back bending of metal films 
experimentally . 24,25 However, the discussion of the exten- 
sion of the experimental techniques from metal films to metal 
nanowires is out of the scope of this paper. For the suggested 
parameter fitting, equation (1) can not be used directly due 
to the lack of Im[k z ], which can not be measured experimen- 
tally. Equation (7) is suitable for the fitting because equation 
(7) is just only valid for Re[k z ] regardless of Im[k z ], This is 
the advantage of our result for the parameter fitting. 

The SPPs back bending at the complex - k dispersion re- 
lation is considered to be induced by the metal Ohmic loss 
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FIG. 2: k : o at the back bending point as functions of (a) the radius r 
and (b) the bulk electron relaxation rate t of the metallic nanowire. 
t = 0.039062 is for the calculation in Figure 2 (a) and r = 0.2 in 
Figure 2 (b). Dielectric constants take the value of 6o = 1 and e„ = 1. 



[ i5 1 i7 1 i8,2i- - 23 i 26 s uc jj conclusion has been drawn in the case 
of planar metal films, but never confirmed in the case of 
metallic nanowires. We show in the following how the metal 
Ohmic loss induces the back bending by using the equation 
(8) for clarity. For a nanowire fabricated with a perfect metal 
r = 0, the denominator of equation (8) has only one term 
of 2(4w 2 - 2). The denominator is real since the frequency 
is real in the complex - k solution. When the frequency to 
approaches the Surface Plasmon frequency co sp — 1 / y2, the 
Re[k z ] behaves as an asymptotic curve to infinity with the de- 
nominator going to be zero. However, when the metal loss t is 
introduced, the denominator in the equation (8) will no longer 
be a real value, but a complex, which can be found in equa- 
tion (7). The absolute value of the complex denominator will 
never be zero no matter what value the real frequency a> is. 
Thus, the Re[k z ] will be finite instead of infinite at the Surface 
Plasmon frequency of a> sp . This is the key point to understand 
how the metal loss changes the asymptotic curve of the SPPs 
to a back bending. 

Last, we consider the general case with Eq and e m intro- 
duced. Equation (12) shows that for a perfect metal with r = 
used in nanowires, the dispersion relation is an asymptotic 
curve with Re[k z ] approaching infinity when the frequency co 

is close to (jj sp = J s ^ . When the metal loss is introduced, 
the dispersion relation then has a back bending near co sp . The 
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Re[kJ 

FIG. 3: back bending position of the Re[k z ] dependent on Eo and e m 
for the metallic nanowire with r = 0.2 and r = 0.039062. Figure 
3 (a) shows the dependence of the back bending position on Eo with 
Eoo = 9. Figure 3 (b) shows the dependence of the back bending 
position on e„ with Eo = 1 . 



dependence of the back bending on the dielectric constants 



8o and 8oo is shown in figure 3. The comparison between the 
analytical solution of equation (12) and the numerical solu- 
tion of equation (1) has been performed, confirming that equa- 
tion (12) can be used to describe the back bending accurately. 
However, for clarity, we only show the analytical results by 
solid lines in the figure. Figure 3(a) shows the dependence of 
the back bending position on Eo with the fixed value e M = 9, 
indicating that changes to be large with the increase of eo- 
Figure 3(b) is for the dependence on Eoo with the fixed value 
8o = 1 , which exhibits that k A) changes to be small with the 
increase of 8oo. The dependence of the back bending position 
on dielectrics shown in figure (3) indicates that equation (12) 
can also be used for the parameter fitting of the dielectric con- 
stants of the medium and the metal. The dependence of the 
complex - to solutions on the dielectric constants is not shown 
due to the tedious reason, which has an asymptotic curve with 
Re[k z ] approaching infinity when the frequency to is close to 

V. SUMMARY 

We have derived approximate analytical solutions to the 
SPPs dispersion relations with n — of cylindrical metal 
nanowires for the relation sections whose frequencies are 
close to the Surface Plasmon frequency in this paper. The 
back bending of the complex - k dispersion relation has been 
confirmed to be originated from the metal loss, which intro- 
duces an imaginary value into the denominator of the expres- 
sion of Re[k z ] and moves the value of Re[k z ] from infinity to 
finity at the Surface Plasmon frequency, resulting in the back 
bending of the dispersion relation. The analytical solutions 
are suggested to be used for the parameter fitting of the phys- 
ical parameters of the nanowires by taking the advantage that 
Im[k z ] can be unknown in the fitting. 
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